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Abstract 



We employ the LSZ reduction formula for Matrix Theory introduced in our 
earlier work to compute the i-pole S-matrix for three form-three form scat- 
tering. The result agrees completely with tree level D = 11 SUGRA. Taken 
together with previous results on graviton-graviton scattering this shows that 
Matrix Theory indeed reproduces the bosonic sector of the D = 11 SUGRA 
action including the Chern-Simons term. Furthermore we provide a detailed 
account of our framework along with the technology to compute any Matrix 
Theory one-loop t-pole scattering amplitude at vanishing p~ exchange. 
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1. Introduction 

It is now commonly believed that eleven dimensional supergravity Q is the low- 
energy effective theory of a more fundamental microscopic theory, known as M- 
theory 0. A non-perturbative definition of M-theory has been conjectured to be 
given in terms of the large N limit of a quantum mechanical supersymmetric U(N) 
Yang-Mills model called Matrix Theory |||. Since the time of this conjecture, many 
computations in various settings have been performed to test the proposal []. Most of 
these works, however, involve the comparison of classical gravity source-probe actions 
with the background field effective action of the super Yang-Mills quantum mechani- 
cal system evaluated on straight line configurations. Clearly, however, a principle test 
of the conjecture || would be to compute the tree level S- matrix of D = 11 SUGRA 
in Matrix Theory. We began this project in where we found, using a formalism 
which enabled the computation of true scattering amplitudes in Matrix Theory, that 
indeed the D = 11 graviton-graviton tree level t-pole S-matrix agrees precisely with 
that obtained from Matrix Theory. We stress that what was computed was the full 
field theoretical amplitude, i.e. some 66 terms depending on physical polarisations 
and momenta, in contrast to previous works yielding phase shifts in semi-classical 
eikonal scattering. 

In this paper, we turn our attention away from the pure gravity sector of the 
theory and consider three form scatteringF] Again, making use of the leading effective 
potential for D-particles at one-loop ||, [7|, §, we find complete agreement between 
the Matrix Theory and D = 11 SUGRA S-matrices (at tree level for the i-pole), an 
amplitude consisting of 103 terms. Together with our previous results on graviton- 
graviton scattering |J this computation confirms that Matrix Theory describes all 
bosonic three-point interactions of the D — 11 supergravity action, including the 
Chern-Simons term / dC A dC A C . In this sense D = 11 SUGRA emerges as the low 
energy effective action of Matrix Theory. 

In addition we give a detailed account of the formalism which could only be 
sketched in our previous letter 0, allowing one to compute any t-pole zero p~ ex- 
change S- matrix element in Matrix Theory at one loop. Finally, as in ||, throughout 
the paper we work in the N = 2 sector of the matrix model, so that we are considering 
the Susskind finite N generalisation [[J of the Matrix Theory conjecture. 

The paper is organised as follows. The main idea of our framework is that S"-matrix 
elements can be constructed from the asymptotic particle states of |10| and involve 
the expectation of the usual Matrix Theory effective potential (including background 

1 Scc Q| for an exhaustive list of references. 

2 The three form contribution to the linearized D = 11 SUGRA potential has been computed in 
Matrix Theory by || . 



1 



fermions) between in and outgoing polarisation states. We refer to this relation be- 
tween Matrix Theory effective actions and ^-matrix elements as the Matrix Theory 
Lehmann, Symanzik and Zimmermann (LSZ) reduction formula and give a detailed 
derivation of this result in section two. Although the one-loop Matrix Theory effec- 
tive action could be derived using quantum mechanical Feynman graphs it is obtained 
more efficiently from the one-loop effective potential of two moving D-particles, com- 
puted using the Green-Schwarz boundary state formalism 0, ||. In section three 
we present a systematic derivation of this result. In section four we combine the re- 
sults of sections two and three to compute the Matrix Theory three form-three form 
scattering amplitude. Furthermore we include a detailed account of the Matrix The- 
ory kinematics involved, along with the algebra of spinor bilinear operators acting on 
polarisation states, an essential ingredient for efficient computations of such ampli- 
tudes. Section five presents the D = 11 SUGRA tree level computation of the three 
form-three form scattering amplitude plus its reduction to the vanishing p~ exchange 
kinematics described by elastic N = 2 Matrix Theory scattering. Finally we give 
our conclusions and in an appendix spell out some conventions and Fierz identities 
needed for our computations. 

2. LSZ for Matrix Theory 

The purpose of this section is to provide the link between the canonical operator based 
Matrix Theory scattering amplitude calculations of ]l(| and the "standard" path 



integral, background field, effective action approach of, e.g., [[12], [13|. This provides 
the technology required for scattering amplitude calculations in Matrix Theory. 

2.1. Asymptotic Particle States in Matrix Theory 

We begin with a short review of the asymptotic state analysis of [I0|. We shall 



primarily study the U(2) Matrix Theory Hamiltonian whose coordinates take values 
in the adjoint representation of U{2) i.e., 

X m = X° m il + X A ia A m = l,...,9 

e = e°ii + e A ia A (2.1) 

where a are the Pauli matrices. Employing a vector notation for the SU(2) part in 
which X m = (X^, X^, X^) and similarly for 9, the Matrix Theory Hamiltonian is 
given byf] 

H = ±P°P° + (|P m ■ P m + \ (X m x X n f + I X m ■ 9 lm x 6) . (2.2) 



3 Note that we are using a real, symmetric representation of the SO(9) Dirac matrices in which 
the charge conjugation matrix C = 1. 
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It must be augmented by the Gauss law constraint 

L = X m xP m - l -9x9 (2.3) 

which annihilates physical states. We wish to study particles widely separated in (say) 
the ninth transverse direction, it is therefore useful (but not necessary), following ||14|| , 
to choose a frame where (say) Xg lies along the z axis 

X] = = Xl . (2.4) 

Denoting the Cartan coordinates x m = X^ and the remainder Y* = (with I = 1, 2 
and a = 1, ... ,8) one then sees that the Hamiltonian (p.2|) takes the form 



H = H CoM (X°) + H v (x m ) + H uo {Yj, J M + H 4 (Y a T , x m , 9 1 , 6 3 ) (2.5) 

Here Hq m = — \ i.dx%) 2 is the U(l) center of mass Hamiltonian. In particular, one 
now observes that Hy = — ^ (d Xm ) 2 xg represents free particle propagation along 
the "Cartan valley" , whereas i^HO describes a system of 16 superharmonic oscillators 
transverse to the Cartan valley with frequency r = (x m x m ) l l 2 depending on the 
distance from the valley origin. Finally constitutes all remaining terms. 

In it was shown that in the limit Xg — > oo, which one interprets as the large 
separation in the nine direction between a pair of asymptotic particles (importantly, 
observe that in a general frame, xg = (Xg ■ Xg) 1 / 2 is gauge invariant, i.e. commutes 
with the constraint), there exists a split of the Hamiltonian into a free and an interact- 
ing part (the latter of which is suppressed in the large xg limit). The free Hamiltonian 
admits eigenstates of the form 

\p l m , H^pl, n 2 ) = \0 B , F ) Xm i e^-P^e^+r 2 )-^ \H l )e^ \H 2 ) e0 -et ■ (2.6) 

Here \^B^F)x m denotes the superharmonic groundstate of i^no with vanishing zero 
point energy. Note that the oscillator states depend explicitly on the Cartan coordi- 
nates x m through their frequency r. In the above, p m are the momenta of the two 
particles and their possible polarisations Ti. 1 are those of the graviton, three-form 
tensor and gravitino represented by states 

\h) = h mn \-) mn , \C)=C mnp \-) mn ^ |^> =<|->-, (2.7) 



whose explicit form was constructed in 1 10|| . Finally, the subscripts 6° ± 8 3 in (|2.6j ) 



indicate from which fermionic variables the polarisation states are built. The state 
( |2.6|) is free in the asymptotic limit, i.e., 



Jim#|pL ^pl, H 2 ) = \ [(p 1 ) 2 + (p 2 ) 2 } \p m ,H\p 2 w H 2 ) (2i 
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for large particle separations Xg. Moreover, the state ( |2.6| ) indeed satisfies the physical 
state condition L\p] n ,'H 1 ]p 2 m ,'H 2 ) = 0. 

2.2. The LSZ Reduction Formula for Matrix Theory 

We now turn to the computation of scattering amplitudes and derive, in particular, 
the formula relating Matrix Theory effective action computations with the scattering 
matrix. For a 1 + 2^4 + 3 process one starts with the S'-matrix element 

S fl = (p\H A ;p 3 ,H 3 \eM-iHT}\p\H\p 2 ,H 2 ) 

= J d 9 X'° Anx' 9 2 d 9 x' d 9 X° Anx 9 2 d 9 x ± e -V-P 8 )m^ e -i(M«).»X' 



pU\ M ^H 3 \ x/ JOB,OF\eM-iHT)\OB,OF)a> m \n x ) \n 2 ) 

gO + g$ eO-gS m 00+03 00-03 



J-9 



e HP -P Jm« TOe HP +P Jro^m (2.9) 



where T is the large time during which the process takes place. The measure factors 
47rxg 2 appear in the integral because of the choice of frame ( |2.4|) . Moreover, since 
we are interested in eikonal kinematics, we have chosen asymptotic states describing 
particles widely separated in the nine direction for both the in and out states. More 
general configurations may be handled via the insertion of a rotation operator in (|2.9p . 

For complete clarity, we note the following. Strictly speaking, one should compute 
Sfi = (out| exp (—iHT) | in) with H being the Hamiltonian ( |2.2j ) in a general frame 
and the asymptotic states |in) and (out| written in a manifestly gauge invariant way 
without fixing the frame Xg = 0, as shown in ||10|| . Now as L\in) = = (out|L and 
[H, L] = we have 

Sfi = (out | exp(-zFT)|in) = (out|IIexp(-iftT)n|m) (2.10) 

where II = (volsu^))" 1 J dflsu(2) exp(iQ ■ L) is the projector onto gauge invariant 
states. Therefore one is able to choose a frame X'g = = Xg at both start and 



endpoints and in this way arrives at ( |2.9| ) in which the variables of the fixed frame ( |2.4j ) 
appear. 

The vacuum to vacuum transition amplitude 1^(0^,0^1 exp(-iHT) |0£,0i?) x . m is 
now the object of interest and may be represented as a path integral with appropriate 
boundary values for the Cartan coordinates 

jT{ Xm ,x' m fi*)-iH CM = ^(o B ,0 F | eM-iHT) \0 B ,0 F ) Xm = 

AM D w YD 16 9 I D 9 x m Unxg 2 ) exp(i / dtL) (2.11) 

Jx m (-T/2)=x m ,e*(-T/2)=e* J-T/2 

Here the Lagrangian L is the usual one ( u pq — H" ) associated with the Hamiltonian 



( p.5|) in the special frame (Xg =0). In particular, as a result of this choice of frame we 
have a measure factor 47rxg 2 at each point in the path which may be exponentiated via 
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ghosts. M denotes some normalisation factor]]. Importantly, note that the boundary 
conditions for the oscillator variables Y£ and 9 1 are zero at start and endpoints since 
we are considering a transition between harmonic oscillator ground states in these 
variables. 

It is also essential to observe that the transition element appearing in ( |2.9| ) de- 
pends on operator valued 9 3 , whereas in fl2.11|) one computes a c-number valued path 
integral. To elevate this result to an operator, as needed in the rest of the compu- 
tation, one makes only errors proportional to the square of the momentum transfer 
q 2 which we anyway neglect in this work since they correspond to contact terms not 
detectable in our /JO-brane computation. To see that the 9 3 boundary conditions 



in the transition amplitude ( 2.11 ) are correct, one can change from the sixteen real 
variables 9 3 to eight complex ones and perform a coherent state analysis similar to 



that of |6|. 

We now make the following observation. Firstly, consider the BRST gauge fixed 
path integral of the ten dimensional super Yang-Mills dimensionally reduced to quan- 
tum mechanics 

e iT{x m ,x m ,e i ) = j\f DAD 9 X m D 6 9 DbDc exp(i / dt L BRST ) 

Jx m (-T/2)=x m ,e3(-T/2)=e3 J-T/2 

(2.12) 

where b and c are 577(2) ghosts, Lbrs(A X, 9, b, c) is the dimensionally reduced super 

Yang-Mills Lagrangian with ghost terms and the gauge field A is the time component 

of the ten dimensional Yang-Mills field. This is the path integral considered in most 

Matrix Theory computations (including the boundary conditions quoted in ( |2.12|) ). 

Then, if one takes the gauge choice A 3 = = Xg and additionally integrates out the 

ghosts and remaining gauge field components A 1 (yielding, respectively, the measure 

factor 4tcx 2 and the frame fixed Lagrangian), one obtains exactly the path integral 

( P- 11| ). The ramifications of this simple observation are clear; one may now obtain 

S'-matrix elements from the effective actions produced by the usual Matrix Theory 

computations found in the literature M . The path integral ( |2.12 ) can be computed by 

expanding about classical trajectories X^ = x^(t) = b m +v m t and constant 9 3 (t) = 9 3 

which satisfy the quoted boundary conditions for impact parameter b m = (x' m + x m )/2 

and velocity v m = (x' m — x m )/T. A quantum mechanical Feynman diagram expansion 

in the gauge of one's choice then leads to the required effective action although we 

found it more efficient to exploit the connection between the Matrix Theory and String 

4 In fact, the observant reader will note that in what follows, we assume that the normalisation 
factor Af behaves as (x' 9 xg) -1 to cancel the measure factors of the initial and final integrations over 
the valley coordinates x m . That this is the case can be argued by SO(9) covariance of the final 
result. We also do not compute the overall normalisation of the path integral. Such technicalities 
should, in principle, be rigorously calculable via a careful construction of the path integral similar 
to that presented in p5[. 
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Theory DO-brane dynamics in order to obtain T(x' m ,x m ,6 3 ) (see section three). Our 
LSZ reduction formula for Matrix Theory, relating the effective action to S-matrix 
elements is therefore 

S fl = f d 9 X'°d 9 x' d 9 X°d 9 x 



exp 



( - i{p A - p 3 ) m x' m ~ i{Pi + P3) m X'° + i(p x - p 2 ) m x m + i(p x + p 2 ) m X^) 
(HA (H 3 \j T ( Xm '<' e3 )- iHc ° MT \Hi) \H 2 ) ■ (2.13) 



Finally, we complete this section by noting that the generalisation of this formula to 
higher N Matrix Theory, i.e., SU(N), N particle elastic scattering with vanishing 
p~ momentum exchange is straightforward. One needs only replace the pairs of 
incoming and outgoing polarisation states in ( |2.13| ) by a set of iV such states. The 
effective action becomes that of an SU(N) Matrix Theory computation depending 
oniV-1 Cartan coordinates and the momentum plane waves become those of iV in 
and outgoing particles. 

3. String Computation of the DO-Brane Effective Potential 

In this section, after a very brief review of the Green-Schwarz boundary state for- 



malism [II]], we give a detailed account of the computation of the one-loop Matrix 
Theory potential (see equation (8) of [[]]), first performed in || ||. We consider here 
the D-particle case, but it is clear from |8| that our result is trivially extendable to 
generic p-branes. 

Dp-brane defects W7\ can be described in String Theory by suitable boundary 



states implementing the usual Neumann- Dirichlet boundary conditions, both in the 
covariant [ITR O, |2T1] as well as the Green-Schwarz formalism HTT1. In the latter frame- 



work, the boundary state describing a single flat D-brane is defined by the BPS 
condition 

Q a + \B)=0, Q d + \B) = 0, (3.1) 

where Q\ , Q\ are suitable linear combinations of the 5*0(8) supercharges Q a , Q a , Q a , Q a . 
The solution for \B) turns out to be 

\B) = exp Y, (-M tJ al n a 3 _ n - iM atl S\S%) \B ) (3.2) 

n>0 ^ U ' 

\B ) being the zero mode part 

\B ) =M ij \m-iMja)\b) (3.3) 
with My, M a a, Maa definite SO (8) matrices [11], |7j, depending on the dimensionality of 



the brane 0. In this gauge, the ± light-cone directions satisfy automatically Dirichlet 

5 In writing M 0( j we have implicitly chosen to work in the IIA theory, relevant for the analysis of 
D-particles. 
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boundary conditions, meaning that they are, effectively, Euclidean branes. One might 
think that defining the boundary state for moving branes by simply boosting the static 
22j would then be problematical, however, as explained in J7|, it is possible 



one 



to overcome this difficulty by identifying one of the 50(8) transverse directions with 
the time direction. Thereafter one deduces the corresponding 50(1,9) expressions 
and performs a double analytic continuation to the final covariant result. 

Given these preliminaries, one may then compute arbitrary one-loop n-point func- 
tions of vertex operators V\, . . . , V n 

A n = f°°dt (B, ^ e -^ta'p+(p-- l d/dx+)y i _ _ _ y\ B ^ /g^X 
JO 

with p-=&£ + osc. the light-cone Hamiltonian(the term id/dx + just implements 
the p~ subtraction needed to obtain the effective Hamiltonian in this gauge). The 
configuration space boundary state \B,x) is given by 

\B,x) = (2^'f- p j -^h el<i ' S \ B ) ® l«) ( 3 - 5 ) 

with (q\q') = vol p+ i (27r) 9 ~ p §( 9 ~ p \q — q') and vol p+ i is the space-time volume spanned 
by the p-brane. Equation (|3.4| ) describes the interaction of D-branes, considered as 
semiclassical heavy objects, with n arbitrary states, described by the vertex operators 
Vi. We are interested in computing the leading one-loop effective action T(v, 9 3 , r) of 
two moving D-particles with relative velocity v at their minimum separation r, that 
is, the usual v 4 /r 7 term plus all other terms, bilinear in fermions, related to it by 
super symmetry. Correspondingly, we will consider the following correlator, encoding 
in particular the abovementioned terms 

V=- [°°dt (B, x = Q\ e -^'p + {P~-id/d*+) e {vQ~+vQ~) e v B |£ ? ^ = ^ (36 ) 
2 jo 

Q-, Q- being the S0(8) supercharges broken by the presence of D-branes and Vb the 
vertex operator that boosts the branes to a relative velocity t>j, given explicity by 

V B = Vi £ da (x^d a X^ + ^S-f U S^j (3.7) 

where the direction 1 will be Wick rotated to give the time direction. The factor 1/2 
in ( |3.6| ) has been introduced in order to normalise the v A term to one. 

A configuration of parallel branes preserves 1/2 of the supercharges; in light-cone 
gauge this implies that among the 16 linearly realised supercharges Sfi, Sq, eight of 
them are left unbroken. Equations ( |3.4j ) and ( |3.6| ) require then the insertion of at 
least eight zero modes (that, due to the constraint (|3.1[), can be always chosen to be 
Sg) in order to get a non-vanishing result. This is precisely the number of zero modes 
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provided by the t> 4 /r 7 term and all its related fermionic terms, which are therefore 
determined by massless string excitations alone 0,0- I n this configuration we can 
then consider only the massless part of ( |3.6| ) which simplifies dramatically. Integrating 
over the cylinder modulus t, we obtain 

V = T Tq J ^^{B \e^ -^e< F) \B ) (3.8) 

where are just the zero mode part of the given supercharges, V B F ^ is the 



fermionic part of the boost operator (|3.7f ), T = v /7r(47r 2 a') 3 / 2 is the tension of 0- 



branes, T is the overall time in which the interaction takes place and q spans the 
directions ±,2,.. .,8. The bosonic part of Vb induces conservation of momentum that 
reads q 1 = q ■ v; we fix in ( |3.8| ) and in the following q 1 = q ■ v = 0, which simply 
means we are computing the effective potential at the timeQ t = 0. It is now conve- 
nient to write the So zero mode trace in terms ofB^ = (S'oTalS'o) / '4. Expanding the 
exponentials we find, following ||, 



n,m: 
n/2 + m=4 



9 m r rl 9 a p 1 ^ 
V = V — / q +l»l-limJUa-..Jn-Un 

^ n\m\J (2vr) 9 q 2 



Vn-ViJ^V + ^) [ai -(^ + ! ^) a " 1 7^ 2 2 ..-7^: 1 1 i: (3.9) 
where t n " As is the usual eight-tensor 

t h...i 8 ^ ^ R hi2 R ^H R ^ R iris (3 10 ) 

whose explicit form can be found, e.g., in the appendix of chapter nine, volume II of 



|23|| . Although the explicit computation for the v A ,v 3 and v 2 cases has been already 
performed in ||, we will report them here for completeness 
f 4 -term: 

§7t»««Vi = g.(S) ✓ (3.ii) 

— * 

where Gd(b) is the propagator for a scalar massless particle in d-dimensions. 
f 3 -term: 



= 2 t; V ^^T^W^ + 2 g„r/7 lm ^ + |^7 lm ^) (3.12) 

where qf c = Y^=2QiQi an d Q ,± = 9° ± g 9 ; notice that q 2 = qf c — q + q~ implying 

that q~ = qf c /q + , modulo contact terms that are vanishing in our configuration of 

6 From now on for simplicity we will omit the overall time of the interaction T and the 0-brane 
tension Tq. 
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separated O-branes. It is trivial to verify that the 5*0(1, 9) expression for the term in 
parenthesis in fl3.12| ) is ipY lrrin il) q n (our conventions for the Dirac matrices are given 



in the appendix), with n = ±,2, ...8, which after analytic continuation (that sends 
also v l — * iv l ) leads to 

V 3 = -2tv 2 v m J 0mn d n G 9 (b) = 2iv 2 v m {6 1 mn 6)d n G 9 {b) (3.13) 

where J^ up = ^T^ vp ^, * is the Majorana-Weyl spinor \1> = m, = {v a ,V & ) and 
m, n — 1, 9. We have written the result also in terms of the 5*0(9) spinor which 
is the most convenient way to connect it to Matrix Theory. It is clear that we can 
from now on restrict our attention to terms with no net power of q + ; as shown in 
( 13.121 ), these terms contain enough information to reconstruct the covariant form of 
our amplitudes.. 
v 2 -term: 



where we defined 



V 2 = —J ^—^*™* ViVjUhMp ( q ,r,) (3.14) 



and q = qrf. It is not difficult to see that 



t liljlmn Pv . v . Uimnp ^ ??) = Jlmq Jln ^ + j™»jn 3 ~^ )qm ^ ^ 



with latin indices labelling the indices i,j,--- = ±,2, ...8, whereas greek indices run 
over all 10 directions. Strictly speaking the equality ( |3.16| ) holds only for m,n — 
2, 8. The cases where m, n — ± are given by the terms in ( |3.9| ) with non- vanishing 
powers of q + , as can be easily checked. This is the sense in which ( |3.16| ) and similar 
identities that will follow should be understood. We can now analytically continue 
the right hand side of ( |3.16| ) leading to 

v 2 = \v 2 {J mQq J\ q + J m ^J\ 3 vv)d m d n G 9 {b) = -2{eY m e){e^ n e)v pVq d m d n G 9 {b) 

(3.17) 

where the 50(9) expression follows after a Fierz identity. 
f 1 -term: 

Vl = ^/(Sj?' 1 *"*"'^-*' 9 ''' (X18) 

This case, as well as the next one, is a little more involved, since we have a new 
contribution 

t Hn "' 36 ViU) h ,,j 6 {q,rj) = 12u 1 i jkkj (q,r]) + 24u ljjkki (q,r])-l/2e lm --- Ja uj jl ,,, j(i (q,r]) (3.19) 
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The first term on the right hand side of ( |3.19| ) is vanishing due to the Fierz identity 
( |A.5| ). The second term is 

L0 ljjkki (q, rj) = Ljinij^mjvip^q^ (3 20) 

Using the relations (|A.8|) reported in the appendix, it is possible to verify that the 
SO(l, 9) expression for the e-term is, up to a numerical factor 



q^ v^ 8 q a qp J^ a J^^J^^ (3.21) 



The expression ( |3.21| ) can be brought into the same form as the right hand side of 
( |3.2U| ) using the identity ( |A.6| ). In order to fix the relative coefficient between the two 
non- vanishing contributions coming from ( ^.19[ ), it is much simpler to consider the 
term proportional to (q + ) 3 , in which case the spinor algebra simplifies considerably. 
This term is proportional to 

- 1/2 e lijl --- j %r]j jlj2 r))(r)j j3j4 r])(r]j j5j6 r)) + 24(r)j lj T])(r]j jk r))(r)j ki r]) (3.22) 

By using the identity (|A.7|) , the e-term in ( |3.22|) becomes 8(7]^ r))(r)^ k r)) (777^*77). 
Putting all the results together we find 



Vi = ^ J 0nl J lfl m J^ p d m d n d p G 9 (b) = -jV t (9Y m e)(6r l 0)(9Y l 0)d m d n d p G 9 (b) 



(3.23) 

where the second identity in ( |3.23|) follows from the first one by 50(9) Fierz identities. 
f°-term: 

7f) r d 9 a 

= ^ / iA- q -rt n --^ 1 ..,M, n) (3-24) 



(2tt) 9 q 2 
where 

t h - is u h .. Aa (q,T]) = 24u ijjkkUi {q, V ) - 1/2 e 1 ' 1 -^...^, rj) (3.25) 
The 50(1, 9) expression for the e-term is 

e Pl ... m ^ 10 q a q^J fll ^ Q J^ fli(3 J^ fl(i " f J^ f18 ^ (3.26) 

whereas 



C0i jjkkui (q, V) = ^■J flum -V n -V P -J p a q Qm q n q P q q (3.27) 



1 

70 

Again, the expression ( |3.26| ) can be cast in the form appearing on the right hand side 
of ( |3.27| ). By looking at the (g + ) 4 -term, similarly to the previous case, it turns out 
that the e-term in ( |3.25| ) gives a contribution equal to Su i jjkkiii( y q,f])- We then obtain 

32 

v = ^J^ m J^ n J/ p J p(T q d m d n d p d q G,(b) 

= ^{e^e^er^e^e^H) d m d n d p d q G 9 {b) (3.28) 
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where once again £0(9) Fierz identities have been used to write the second identity 



in O). 

Collecting all terms we obtain the final result for the leading one-loop potential 
of two DO-branes 01 



v 4 + 2iv 2 v m (e 1 mn e) d n - 2 Vp v q (dY m 0)(0i qn e) d m d n 

An 

--v^e^e^e^e^e) d m d n d p 

y 

+^(6 1 ml 6)(6 1 nl 6)(6 1 pk 6)(6 1 'i k 6) d m d n d p d q ] G 9 (b) 



(3.29) 



The first, second, third and last terms of (|3.29| ) were calculated in a super Yang-Mills 

respectively. The supersymmetry parameter 9 should 



context in pi,|2l,p5| and| 
be identified with the spinor 6 3 /2 introduced in the previous section and represents 
the fermionic background in Matrix Theory. 

Before concluding the present section we want to make some comments about the 
origin of the e-terms in (|3.19j) , fl3".25| ). By performing an analysis of 1-point functions 
of massless closed string states on a disc with supercharges inserted on its boundary, it 
is straightforward to derive which fields exchanged between the branes are responsible 
for the interactions described above |27|, 0, §]■ I n this way, as one might expect, all the 
interactions, except those coming from the e-terms, are due to exchange of dilatons, 
gravitons and Ramond-Ramond (RR) vector gauge fields. On the other hand, the 
e-terms arise from an interesting coupling between dual RR gauge potentials, very 
similar to that analysed in |28|, ^9[ [K|.In particular the e-term coming from the part 
of the potential linear in the velocity is due to exchange of a RR one form Am and its 
dual seven-form Am- The insertion of six supercharges corresponding to D-particles 
on the boundary of the disc indeed induces a non-minimal coupling to the seven-form. 
Schematically, in light-cone gauge it reads 



(B\Q 6 \A ( 7)) 



■Ja. 



(3.30) 



where the direction 1 satisfies the Neumann boundary condition on the disc and will 
be identified with the time direction. The coupling (|3.30| ) produces then an interaction 



(B\Q 6 \A (7) )(A {7) \A (1) )(A (1) \V B \B) ~ ± Vl u n . 



■J a 



Hji-je 



(3.31) 
: 51(7). 



where we used the relation between dual transverse gauge potentials A(i) 
Notice the relationship of this interaction to the corresponding one analysed in |30 
for the D0-D6 brane system. 

The e-term associated to the static potential can be treated similarly. Again, the RR 
seven-form Am has a non-vanishing 1-point function when eight supercharges are 
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inserted, that includes a term 



(B\Q s \A (7) )^u; ln ..,MAh-M (3-32) 

leading to an interaction 

(B\Q S \A [7) )(A (7) \A {1) )(A {1) \B) ~ ^ J ,.. i8 (g)e^'^ (3.33) 

equal to the e-term appearing in the static potential. In this case the correspondence 
with the analogous static RR potential for the D0-D8 system, found in f29R , is more 



subtle. Indeed, it has been shown in that non-physical polarisations of the RR 



nine-form, identified with some RR one-form polarisations, are responsible for the 
RR attraction between a DO and a D8 brane. Although these effects are clearly not 
visible in a transverse physical gauge, it can be seen in a covariant formalism that 
the RR nine-form has a non-vanishing 1-point function with DO-branes, when eight 
supercharges are inserted. The considerations above are of course only schematic but 
the key point is to highlight the presence of such interesting interactions which can 
be analysed in more detail along the lines of [25|, ^9, ^ . 



4. Three Form Scattering in Matrix Theory 

In this section we present the results for three form scattering in Matrix Theory. 
The section will be divided into three parts. In the first, we spell out carefully 
the kinematics of the scattering amplitude under consideration. In the second we 
develop the algebra of bilinears built from the 5*0(9) fermionic operators 9 3 acting on 
polarisation states. This latter development allows one to perform the Matrix Theory 
computation of S- matrix elements, once one is given the potential as in Q3.29j ), with 
comparable (or even improved) efficiency to that when employing tree level Feynman 
diagrams in D = 11 supergravity. In the third section we state our result. 

4.1. Kinematics 

The starting point is our Matrix Theory LSZ formula (|2.13|) which yields the S- 

matrix for the 1 + 2 — > 4 + 3 scattering of particles with momenta p x ml p 2 ml p^, 

and p^, respectively. To begin with, the free U(l) center of mass sector of the theory 

decouples and yields an overall factor, {^In) 9 5 9 (P' m — P m ) exp(—iP m P m T/2), expressing 

conservation of total momentum. Here P m = p^+p^ and P' m = Pm + p^ are the total 

in and outgoing momenta, respectively, and the exponential is the standard factor 

obtained in time independent perturbation theory^]. 

7 In what follows, we shall disregard these kinematical prefactors in comparing to the SUGRA 
Feynman graph result since they only express the usual relation between time independent and time 
dependent perturbation theory. 
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A loopwise expansion of the Matrix Theory effective action yields 

T{x' m , x m , d 3 ) = T(b m , v m , 6 3 ) = v m v m T/2 + (4.1) 

with v m = (x' m — x m )/T and b m = [x' m + x m )/2, so that the S matrix now reads (up 
to an overall normalisation) 

S fi = J d 9 x'd 9 xexp(-tw m x' m +tu m x m +tv m v m T/2){H 3 \ {H 4 \ ^ (loops) M^) ^ ^ 

(4.2) 

where we denoted the in and outgoing relative momenta u m = (p m — p m )/2 and 
w m = {.Vm ~Pm)l^i respectively. However, changing variables d 9 x'd 9 x — > d 9 (Tv)d 9 b, 
the integral over Tv m may be performed, for large T, by stationary phase which yields 

S fi = e-^ u+ ^^ 2 J d 9 be~ i « mbm (H 3 \ (ft 4 | ^°°^(v m =(u +w)m /2,b m ,e3) ^ ^ (43) 

where r( loops ) is to be evaluated at v m = (u + w) m /2, which we take, henceforth, as 
the definition of v m . Moreover q m denotes the momentum transfer q m = w m — u m . 
For clarity, let us give the relation between the various momenta introduced above 

Vm P'm Pm (Zm/^ Vm P'm V"m Pm Qm/^ 1 

Vm Pm ^m Pm ~\~ V m ~\~ Q m /2 Vm Pm ^m Pm Qm/^ 1 

Energy-momentum conservation fixes, in particular, the important relation q m v m = 0. 

Finally, we make one last manipulation. The Matrix Theory effective action is 
an expression of the form r( lo °P s ) = dt y(ioo P s) = ^ + Vm t),x m = v m , 6 3 } for 
some potential y ( 1oo p s ) . However, we would really like to consider a time independent 
potential depending on b m and v m and 9 3 acting on states. This is achieved by 
expanding the exponential inside the polarisation expectation value in (fOf) and then 
interchanging the dt and d 9 b integrations. Then shifting b m — > b m + v m t and using 
QmV m — 0, the dt integral yields only an overall factor T (the same factor T as 
appearing in front of equation ( [3.8|) in our string computation). Finally, we have the 
desired expression for the one-loop Matrix Theory S'-matrix 

S^ = iTexp(-iv 2 ) Jd 9 be~ iqmbm (n 3 \ (H 4 \ V (1) (6 m , v m , 6 3 )\H l ) \H 2 ) . (4.4) 

The loopwise expansion of the effective action is valid for large impact parameters 
b m and hence small momentum transfer q m . In section five we will show that this is 
precisely the limit dominated by t-channel tree level physics in the SUGRA compu- 
tation. 

4.2. Algebra of Bilinears 

All that remains to complete our Matrix Theory computation is to perform a 
Fourier transform with respect to b m of ([4.4j) and insert the one-loop Matrix Theory 
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potential ( |3.29| ) into the polarisation inner products. The result is an amplitude which 
may be directly compared with the result of the SUGRA tree level Feynman diagram 
calculation of section five. 

In our previous work, when computing the graviton-graviton scattering amplitude 
we used the explicit representation for the polarisation states in terms of 50(7) £g> £7(1) 
covariant complex spinors (see |]IU|| ). Then rewriting the Matrix Theory potential 
also in terms of complex spinors, we were able to compute its polarisation expecta- 
tion value. Unfortunately, if one considers now the three form polarisation states, 
this naive method becomes rather quickly unwieldy. Instead, we now present an 
50(9) covariant algebra of spinor bilinears acting on polarisation states with which 
amplitudes of complicated Matrix Theory potentials can be efficiently computed. 

The Matrix Theory potential (|3.29|) is written solely in terms of 50(9) rotation 
generators 1 9 3, y mn 9 3 where the operators 9 3 a satisfy the anticommutation relations 
{0a,9p} — da/3- However, the states 

_ 'iyl>4 I \M \n-/2,3\ _ 1/2,3 i \M (a r\ 

\ n I — n M I - /eo+es , \rL ) — rt M \— )qo_qs (4. 0) 

(where the generalised index Ai denotes M. = {inn; mnp;ma}, corresponding to 
graviton, three-form and gravitino polarisations, respectively) are built from states 
depending either on the sum or difference of the centre of mass ({9^,9®} = 5 a p) 
and Cartan fermionic coordinates which we denote as 9 l = (9° + 9 3 ) / v2 and 6 2 = 
(9° - 9 3 )/V2. Of course, the operators \9 1 1 mn 9 1 and \9 2 1 mn 9 2 act simply as 50(9) 
rotation generators on the states |— )q% and |— )$2 , respectively, under which they 
transform in the usual fashion. Therefore, if we could write the potential only in 
terms of these operators the computation would be completely trivial. Yet there are 
cross terms since 

8 ' 4 ' 4 ' 2 ' K J 

Let us now concentrate on amplitudes involving bose particles only (the generalisation 
to the fermi case is simple but not needed here). Clearly, the difficult terms are those 
involving the bilinear 9 lr ) mn 9 2 . However, since between bose states the expectation 
of an odd number of 9 operators vanishes (and the 1 and 2 sectors are independent), 
only an even number of such operators can occur. But products 9 lr j mn 9 2 9 1 Y S 9 2 
may always be rewritten in terms of the 50(9) generators 1 ' 2 ^ mn 9 1 ' 2 or three index 
bilinears 6 |1 ' 2 7 rmip 6 |1 ' 2 via a Fierz rearrangement, of which only the latter cause any 
difficulty. 

Clearly, what is needed then is the action of three index operators 9^f mnp 9 (drop- 
ping the labels 1,2) on states |— ) M . The result is easily obtained by making the most 



general ansatz and fixing the coefficients via the explicit representation given in [10 
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we find 



T mnP0|_)t« = -i2AV3(5 mt \-) npu - I 5 tu \-) rnn P) , (4.7) 

g^mnpQ^_yuv _ i2A^/3 5 mt 5 nU \ -) pV + - e mn P tuvwyz^_yyz (4 8) 

On the right hand side of ( [4.7] ) and ( |4.8| ) one must (anti) symmetrize (with unit weight) 
over all indices according to the symmetry properties of the left hand sides of these 
equations. 

4.3. Results 

Given the algebra ( |4.7| ) and ( |4.8| ) along with the potential ( |3.29| ), only a moderate 
amount of computer algebra |31|] is now required to obtain the Matrix Theory one- 
loop three form-three form eikonal scattering amplitude. Our result consists of 103 
terms and is given by (normalising the v 4 term to unity) 

A = ^ {I^Cm^ + Sv 2 (C u (v,q)C 23 - C u (q,v)C 23 

+v\\ I C 12 (q, q)C 3A + | § C 12 C 34 (g, q) - § C 14 (g, q)C 23 - § C 13 {q, q)C 24 
+| 1 Cu{q, mi)C 23 (g, mi) + || C u {m u q)C 23 {m u q) - § C 14 (q, mi)C 23 (mi, q) 
+ | 1 C 13 (q, m 1 )C 24 (g, m x ) + | f Ci 3 (m 1; g)C 2 4(mi, g) - | Ci 2 (g, m 1 )C , 34 (g, mi) 
+ | | Ci 2 (g, mi, q, m 2 )C u (mi,m 2 ) + | § (7i 2 (mi, m 2 )C 34 (g, mi, g, m 2 ) 
-| Ci3(g, mi, g, m 2 )C 24 (mi, m 2 ) 



-| | Ci 2 (u, w)C 34 (g, g) + | C 12 (v, q)C u (v, q) + | § Ci 4 (u, g)C 23 (w, g) 
-| C 14 (v, g)C 23 (g, u) + || Ci 4 (g, v)C 23 (g, u) + | Ci 3 (v, v)C 24 (q, g) 
-| | Ci 3 (u, q)C 2A (v , g) - 1 1 Ci 3 (g, ^)C 2 4(g, «) — 1 1 C 34 (w, v)Ci 2 (g, g) 
+| 3Ci 2 («, g, g)C 34 + 1 3C 34 (w, g, v, g)Ci 2 
-6C u (v, g, f , g)C 23 - 3C 13 (v, g, u , g)C 24 

-| 9Ci 2 (mi, m 2 )C 3 (u, g, mi)C 4 (?;, g, m 2 ) - | 9C 34 (mi, m 2 )C 1 (v, g, mi)C 2 (u, g, m 2 ) 

+9C 24 (mi, m 2 )Cx{v , g, mi)C 3 (i;, g, m 2 ) 

-9Ci 2 (f , g, mi)C 34 (g, m x ) + 9Ci 2 (u, g, g, mi)C 34 (u, m x ) 

+9Ci 3 (w, g, v, mi)C 24 (g, mi) - 9Ci 3 (f , g, g, mi)C 24 (f , m x ) 

-9Ci 3 (mi, v)C 2A (q, mi, u, g) + 9Ci 3 (mi, q)C 24 (v, mi, v, q) 

+9Ci 2 (mi, v)C 34 (q, m h v, q) - 9Ci 2 (mi, q)C u (v, m 1 , v, q) 

+ | 18Ci 2 (t> , mi, v, m 2 )C 34 (g, mi, g, m 2 ) + 1 18C 34 (v, mi, v, m 2 )C 12 (q, mi, g, m 2 ) 

-18Ci 2 (f , mi, g, m 2 )C 34 (f , m x , g, m 2 ) 

+| Ci 2 (v, g)C 3 4(g, g) - | C u (v, g)C 23 (g, g) + f Ci 4 (g, v)C 23 {q, g) 
+! CisK g)<^24(g, g) - | C 13 (q, v)C M (q, g) - f C 34 (t>, g)Ci 2 (g, g) 
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+| C 12 {v, q, g, mi)C 34 (g, m x ) - § C u (v, q, q, m 1 )C 23 {q, »i) 

+ | C u (v, q, q, mi)C 23 (mi, g) + § Ci 3 (u, g, g, mi)C 24 (g, m x ) 

-§ C 24 (g, mi, d, g)Ci 3 (mi, g) - | Ci 4 (g, mi, g)C 23 (g, mi) 

+ | Ci 4 (g, mi, f , g)C 23 (mi, g) - § C 34 (g, m x , g)Ci 2 (mi, g) 

-9C 12 (f , mi, g, m 2 )C 34 (g, m 1 , g, m 2 ) + 9C u (v, m 1 , g, m 2 )C 12 (q, m 1 , g, m 2 ) 

-| I Ci2(g, g)C 34 (g, g) - ~ 1 C u (g, g)C 23 (g, g) - ~ § Ci 3 (g, g)C 24 (g, g) 

+| | Ci2(g, mi, g, m 2 )C 34 (g, mi, g, m 2 ) + \ § Ci 2 (g, mi, g, m 2 )C 34 (g, m 2 , g, m x ) 

+ | | Ci3(g, mi, g, m 2 )C 24 (g, m 2 , g, mi) 



+ 



d <— C 2 , C 3 <— C 4 } (4.9) 



where we have introduced the notation in which an n-index tensor T mi ,.. mn written 
as a function of a vector g m denotes T(g, m 2 , . . . , m n ) = T mim2 _ mn q mi . Moreover the 
tensors (z, j = 1, . . . , 4) denote the contraction of polarisation tensors, for example 
C 14 (g, v) = C^ im2m ' A CT 3m2ni q mi v Ul (where the indices are contracted between the two 
tensors in the order indicated, another example is Ci 2 = c™ im2m3 C% lsm2mi ). In the 
next section we shall see that this result^ yields perfect agreement with SUGRA 
Feynman graphs. 

5. The Supergravity Computation. 

The bosonic sector of 11-dimensional supergravity [p]] is given by 

£ = — 2^ V-g R- IV~9 \f~mnpq) 2 

~ 12^ £Ml "' Ml1 RmiM 2 m 3 m 4 Fm 5 m 6 m 7 m s Cm 9 M 10 Mu (5.1) 

where Fmnpq = ^9[mCnpq] and g = det gMN- Perturbative quantum gravity may be 
studied by considering small fluctuations h MN from the flat metric rj MN 

gMN = Vmn + K h MN (5.2) 

where k is the 11-dimensional gravitational coupling constant. From now on we raise 
and lower indices with the flat metric t]mn- Propagators are obtained in the usual 
fashion. For the graviton we employ the harmonic (de Donder) gauge <9jv/i jv m — 
(l/2)dMh N n = with propagator (in d dimensions) 

— 4i / 1 \ 

(h MN (k)h PQ (-k)) = -^-[r] iM \ P] T] N) Q - j—^VmnVpq) • (5-3) 



8 We note, in passing, that in obtaining (4.E), products of as many as four nine dimensional Levi- 
Civita symbols were encountered and expanded in Kronecker deltas. For convenience, we have also 
put an explicit i in front of terms mapped to themselves under the replacement [ 1 < — > 2 , 3 < — > 
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For the antisymmetric tensor, the gauge fixing function 3mC m np m the weighted 
gauge £fi x = — \{d M Cmnp) 2 yields the Feynman propagator 



(C MlM2M3 (k) C N ^(-k)) = s$ h sZ Cly (5-4) 

The relevant vertices are easily read off from (BTTI) . In particular, note that Hmn 
couples in the usual way to the three form stress-energy tensor. 

At tree level, the only graphs contributing to four point, three form scattering are 
the single graviton and three form exchange diagrams. These are easily computed 
and in the t = —2p x M p^ channel, which, as we shall see dominates eikonal physics 
(the s and u channels follow anyway by Bose symmetry) one finds 

^84 = ^l{±(eFW)»» p {eF*F*) MJfP 



-32 



(5.5) 



The factor R in (|5.5f) represents the radius of the compactified light-like circle, 
F 1 MNP q = P\m^npq] (i = 1) ■ ■ ■ >4) is the curl of the eleven dimensional polarisation 
tensor C l MNP and we denote F l ■ F j = F i MNOP F^ NOP , (F l ■ F j ) AB = F\ MNP F 3 B MNP 
and ( e F i Fi) MNP = e MNPMx - M * Nl - N *Fj il ^F% v .. N4 . Of course, one must put D = 
11 in this formula. The momenta and polarisations satisfy the mass shell and gauge 
conditions p^pf 1 = and Pi I C t MNP = 0, respectively. 

Since we are not aware of any other occurrences of the result (|5.5|) in the literature, 
for added certainty we considered the dimensional reduction of the amplitude ( |5.5| ) 
to ten dimensions setting p l n = 0. In this case, the three form gauge field, that 
transforms as an 84 with respect to the little group 5*0(9), splits into 56 and 28 
representations of SO(8), that is the little group in ten dimensions. Then, we checked 
the validity of ( |5.5| ) by computing the four-point function of the Ramond-Ramond 
three form (56) and of the antisymmetric tensor field (28) in type IIA string theory. 
By taking then the low-energy a' — > limit of the IIA amplitudes we obtained 
precisely the dimensional reduction of (5J5) to ten dimensions. 



5.1. Light Cone Supergravity Amplitudes. 

In order to make a comparison with the Matrix Theory results we need to rewrite 
our supergravity t channel amplitude in terms of physical transverse nine dimensional 
degrees of freedom. Since we are considering the N = 2 Discrete Light Cone Quanti- 
sation (DLCQ) formulation of the theory, we work in light cone coordinates and 
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specialise to the case of vanishing p momentum exchange. Define, therefore f\ 

„10 I n o 

P =Pt= -j= ( 5 - 6 ) 

so that on-shell momenta satisfy = p Pm — 2p + p~ + p m p m where nine dimen- 
sional indices m,n, . . . are contracted with a Kronecker delta. The gauge condition 
P M Cmnp = may then be solved in terms of 84 physical polarisations C mnp via 

C*H — m C-\-mn > C— mn Pr^ rmn • (^•'^) 

P 

where we have used the residual gauge freedom to set C +mn = 0. From now on, 
we measure momenta in units of the compactified radius R so that p~ = 1 and 

The N = 2 DLCQ Matrix Theory describes 1 + 2 — > 4 + 3 scattering in the case 
of vanishing p~ momentum exchange so we may therefore write the incoming and 
outgoing momenta as 

Pm = ("I ( v m - <?m/2) 2 , 1 , v m - q m /2) p 2 M = (-| (v m - q rn /2) 2 , 1 , -v m + q m /2) 
Pm = (-| ( v m + <?m/2) 2 , 1 , v m + q m /2) p 3 M = {-\ (v m + q m /2) 2 , 1 , -v m - q m /2) 

(5.8) 

where we have set the nine dimensional centre of mass momentum to zero by trans- 
verse Galilean invariance. Importantly, note that conservation of p + momentum im- 
plies v m q m = 0. The Mandelstam variables in this parametrisation read 

t = q 2 m = -2p l M pf , s = Av 2 m +q 2 m = -2p x M p™ u = Av 2 m = -2p l M p^ = s-t . (5.9) 

Eikonal scattering, for which the scattering angle Cos _1 ([t> 2 — q 2 /A\/[v 2 + q 2 /4]) ~ 
q 2 /v 2 is small takes place for small q 2 = t. Therefore we must study the ampli- 



tude ( |5.5| ) in the small t limit which is dominated by the t-pole. In fact, we argued 
above (see section three) that our DO-brane computation was not reliable for contact 
terms proportional to q 2 , and therefore in this work we only study those terms which 
do not cancel the t-pole (one may therefore disregard all terms with an s or u pole). 
We will further discuss this restriction in the conclusion. 

It only remains now to state our main result, namely that substituting equa- 
tions (|5.7|) , (|5.8| ) and ( |5.9|) into the amplitude (|5.5| ) and neglecting terms cancelling the 
t-pole, one reproduces precisely the 103 terms of the Matrix Theory amplitude 



^Notice the conventions used here are slightly different from those used in section three. 
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6. Conclusions 

In this paper we have computed a four-point scattering amplitude and shown, 
analogously to the graviton case analysed in ||, that Matrix Theory also reproduces 
correctly the tensorial structures of the D = 11 SUGRA three form couplings. Both 
the present and our previous matrix computation have been performed by considering 
only the leading one-loop DO-brane effective potential that is protected by supersym- 



metry from quantum corrections |32|]. Therefore the agreement found might have been 
expected as is the case for most of the one-loop phase shift calculations appearing 
in the literature. Furthermore, in principle, it should be possible to fix the structure 



and coefficients of equation ( |3.29| ) by supersymmetry alone |32|]. However, without 



a formalism in which amplitudes can be derived from the Matrix potential ( |3.29| ), it 
would have remained unclear how such a potential could be compared with the ten- 
sorial structure of the SUGRA amplitude. Clearly our work resolves this issue and 
renders the precise relationship between supersymmetry in each of these models more 
transparent. We also remark that our formalism might provide a route to establishing 
eleven dimensional Lorentz invariance of Matrix Theory. 

Another issue deserving comment is the choice of background about which one 
perturbatively expands D — 11 SUGRA in order to compare with Matrix Theory. 
For the case of finite (and actually small) N, it is clearly natural to expand about a 
flat background and ignore the geometry induced by the D-particles themselves (since 
in that case these states appear as fundamental Kaluza-Klein states, i.e., excitations 
of the flat vacuum), in contrast to the large N case when they can be represented as 
classical sources, modifying then the background geometry. 

Throughout this paper we have considered t channel amplitudes only. Since we 
consider quantum asymptotic states which can describe in and outgoing particles 
scattering at any angle, this restriction can in principle be relaxed, although it is not 
completely clear how to compute the transition element fl2.9| ) (or equivalently the 
corresponding path-integral ( |2.12| )) in this case Furthermore, our current work is 
also limited to the case of vanishing p~ momentum exchange; although of considerable 
interest, interactions involving p~ exchange are clearly not visible in our perturbative 
Matrix Theory framework. 

Finally, an important line of development could be to use our technique to analyse 
higher order Matrix Theory scattering amplitudes that would correspond to, say, one- 
loop effects in supergravity. This kind of comparison will put Matrix Theory to a 
much more stringent test and will be crucial in trying to understand the range of 
validity of the theory itself. 

10 In the context of classical gravity source-probe approach, recoil effects have been recently taken 
into account in |33|. 
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A. Appendix 

We employ the following Dirac matrix conventions, the 32 x 32 50(1, 9) matrices are 

r° = ( 7 ° J ^), r=( 7 ^)), , = i,... ,9 (A.i) 

\-i(i6) J J 

where the 50(9) Dirac matrices are chosen to be real and symmetric 

h) o \ , _( o 7 ; 8) 



The 50(8) Dirac matrices rfg) are the same as those appearing in volume one of [23|. 
The charge conjugation matrix C is identified with T so that for a Majorana spinor 
iff = \1/ T r°. In section three use has been made of the following Fierz identities 

t fiy 7 * = (a.3) 

which in 50(9) and 50(8) language read 

(^i)0)(*7?ii)*) - miumiU ) = ( A - 4 ) 

(v&(v7^v) - 2(vi(sMvi l (s)V) - (vilivmfvv) - MMvi&v) = (A.5) 



with iff = (qY 9 = (^ a ,^ a ) and the common notatioiry 11 '" 1 " = (l/n\) 7* 1 ../y tn ztperm. 
valid for all gamma matrices. 

The expressions (|3.21|) , (|3T26|) reduce to the form appearing in the right hand sides 
of equations ( |3.20| ) and ( 13.271) , using the relation 



q a q p if>T^ 2a if> if!T^ 4l5 if ~ q a q p ^r a ^ $r wM2W«^ t + _ _ _ ( A6 ) 



where dots stand for terms proportional to 8^, 8 ft, 8ft that vanish in ( |3.21| ),( |3.26| ) and 
contact terms proportional to 8% have been neglected. In order to write our potential 
in the form shown in (gjgg) , the 50(9) Fierz identity (A.5) of f2f| as well as flAH ) 
are needed. We used moreover the Fierz identity 

( V1 ^ V )(m kl] ?j) = -le^^irj^ifh^f}) (A.7) 
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as well as those involving the e tensor 

(A) ■ q h (VJ i2 fl)(vi i3Um v) ~ (A) ■ q H [( V J i2m v)(vi i3ii V) ~ (Vl^vKvf 2 ^ 

(A) • qk (vY li2k v)(m l ' iUm v) ~ (A) • g l4 [(vi ili2 v){vi mi3 v) + (vi mi3 v)(vi lll2 v)}; 

(A) ■ qil (vi m v)(vi l2l3l4 v) ~ (A) • ftl [(ttt^X^t* 3 **^) - (^NXw^)]; 

(B) • (^^VKVI^V) ~ (B) • (r ?7 n ^)(r/7 i2 ^ 4 ' m ^); (A.8) 

«i ft, g m g« [(^7T iw ^)(w i2 ' 3i4 ' s ^) + 2(^ 7 ia ^) (777^^77)]; 

e n --- 48 gng m gng fc (^7 i2m r/)(^ 3n ^)[(^7 fe ^)(^7 iJ4 --- i8 ^) + (^^(^T* 4 "** 7 *^] 
where 

e iih...i 6 Vi q m q n (r]Y 5i6n rj) 
e n '" tB q m q n q p (r] r y 1 ' 51 ' 6n fj)(r]Y r%8P! n) 

e h ~ ia qh q P q n q m (fjf si6 fj)(r]Y 7isp fj) 

valid, respectively, for the linear in v and static term in the potential. As explained 
in section three, in order to fix the relative factors between the two contributions 
appearing in the linear in v and static effective potential, there is no need to know 
the coefficients in ( |A.6| ) and ( A.8|) , as it is by far more convenient to use only ( |A.7 ). 



(A) 



(B) 
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